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ABSTRACT 


In this paper, the maximum downwind distance Xmax is calculated which depends on the dispersion parameters o, and o, by using 
two methods as follows: the power law function method, and standard method at three effective heights at 5, 100 and 250m 
respectively. The maximum ground level concentration of air pollution that based on the Gaussian plume model can be calculated 
by using the values of the estimated and actual Xmax in the equation of Gaussian air pollution model for different stability classes. The 
estimated and actual Xmax are summarized in power law and standard method. Also the percentage of the error between estimated 
and actual maximum downwind distance are in the range of 0 to 492 and 0.005 to 14.2 respectively. Comparison between the results 
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of estimated and actual Xmax in case of power law and standard method with the previous work used Brigg’s method (Ronbanchob 
Apiratikul, 2015) are done. 


Keywords: Gaussian model; maximum downwind distance; maximum concentration; dispersion parameters. 


1. INTRODUCTION 


The concentration of pollutant resulting from a source or group of sources is the main factor of air dispersion modeling under 
various meteorological conditions. These models are useful for studying the transmission of pollutants into the air (Essa et al., 2015). 
The concentration of pollutant is a function of number of variables, including the emission rate, the distance of the receptors from 
the source, and the atmospheric conditions. The atmospheric air quality dispersion model was usually used to estimate the 
reduction that occurred during the transport of pollutant from an industrial source (wang et al., 2006). 

The most common model for studying the air dispersion process was used to improve air quality based on Gaussian plume/ puff 
formula. The Gaussian model was derived by many authors such as Sutton (1953), Turner (1970) and Smith (1973) was contained 
many parameters on the basis of empirical correlationa, and o, as function of distance were originally developed by pasquill and 
modified by Gifford. The lateral and vertical dispersion parameter, respectively o, and o, represented the turbulent parameterization 
in this approach. They contained the physical ingredients that described the dispersion process (Abdel Wahab et al., 2013). 

In this paper, the maximum concentration at ground level in center-line is predicted with the maximum downwind distance by 
using the dispersion parameters by two different methods. The Power-law function was the most commonly used method in which 
plume dispersion coefficients was expressed in terms of downwind distance and atmospheric stability (Sharan et al., 1995). Also the 
standard method which was based on a single atmospheric stability determined by vertical temperature gradient also the analytical 
expressions based on (P- G) curves (Essa, et al., 2005). A Comparison between the results of estimated and actual Xmax in the case of 
the power law and standard method with the previous work used Brigg’s method (Ronbanchob Apiratikul, 2015) were done. 


2. MATHEMATICAL MODEL 


The equation of Gaussian air pollution model can be expressed as equation (1) (IAEA safety guide, 1983): 


Ip -2) ics / cy 


e az! +e 2laz (1) 


c(x, y,z,H) = 


2MdyozU 


where C(x, y, z, H) is the concentration in the downwind distance, crosswind and vertical distance at stack height H, dp is the 


emission rate at point source "p",ay ,a, are crosswind and vertical dispersion parameters respectively, u is the wind speed (m/s). 


Equation (1) can be reduced to simple term at the ground level in the plume centerline as follows: 


c(x, 0,0, H) = —” les ; 


TOyozU 


where the variables oy , 


be derived from the principle of rough estimation by solving equation o, = = (Kenneth and Cecil, 1972) using estimated and actual 


a, depends on the atmospheric stability classes. The maximum downwind distance and concentration can 


methods, so we used two different methods such as power law and standard methods as follows: 


1. Power-law functions. 
The variables a, and a, can be calculated as follows (Essa et al. 2006): 


Tyaax? (4) 
Tract (5) 


where a, b, c, d are constants values depending on the stability classes (Table 1) (Henderson, 1989). 
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Table 1 Constants for calculating lateral ( [ o] _y) and vertical dispersion parameter ( [ o] _z 


Oy Oz, 

Stability a B c D 

Very unstable 0.40 0.91 0.41 0.91 
Unstable 0.36 0.86 0.33 0.86 
Neutral 0.32 0.78 0.22 0.78 
Stable 0.31 0.71 0.06 0.71 


2. Standard method 


This method is based on a single atmospheric stability which determined by the vertical temperature gradient, = (Green et al., 1980) 


and the dispersion parameters have the form: 


oy. _(5)0, =a (6) 


where r, s, a, p and q are constants depends on the atmospheric stability (Table 2) (Sharan et al. 1995). 


Table 2 A correspondence between, dg , dispersion parameters and Pasquill stability classes 


Pasquill classes A B C D E 

a (om) < -1.9-1.9t0-17 -17to-1.5  -1.5to-0.5 -0.5to1.5 
69252015105 

a (km) 0.927 0370 0.283 0.707 1.07 

s (m/km) 102.0 96.2 72.2 475 33.5 

q -1918  -0.101 0.102 0.465 0.624 

r (m/km) 250 202 134 78.7 56.6 

p 0.189 0.162 0.134 0.135 0.137 


Now, one estimates the values of crosswind and vertical standard parameters in power law and standard methods as follows: 


Table 3 power law and standard method for prediction of o, and a, 


Stability power law standard method 


Class oy(m) o,(m)oy(m)o,(m) 


0.25x 0.102x 
A(very unstable) 0.4x°.91 0.41x°9! aa 


140927 10.927 


0.202x 0.0962 
B(moderately unstable) 0.4x°.9! 0.41x°9! eee TT 


1+y370 10370 


C(slightly unstable) 0.36x°.86 0.33x°86 lial 


1+ 
0.283 0.283 


0.0787x 0.0475x 


D(neutral) 0.32x°78 0.22x°78 7 DRS | O46 
a 0.707 a 0.707 
0.0566x 0.0335 
E(slightly stable) 0.31x°71 0.06x°71 unr 7 ware 
1.07 1.07 
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3. RESULTS AND DISCUSSION 


The actual Xmax of the Gaussian model is determined using programs Mathmatica5 that finds the relation between the maximum 


ground level concentration and downwind distance through three effective heights at 5, 100 and 250 m were evaluated using qp and 
u at 3 g/s and 3 m/s respectively. The percentage of relative error was calculated by 


Relative error= actualxmax—estimated xXmax 


actual Xmax 


The calculation of Xmax which depends on the dispersion parameters can be estimated from Table 4 


Table 4 calculation of Xmax 


XmaxXmax 
Stability class | power law standard method 
— 
H 0.91 H 
A 
0.412. 0.144+2.07H 
1 
B ( H \oo1 H 
0.41v2. 0.136+0.27H 
4 
H \o8e H 
Cc 
0.33V2. 0.102—0.360H 
1 
D (ha) oat 
0.22v2. 0.067-0.657H 
1 
E ( H ae H 
0.06V2. 0.047—0.58H 


The prediction of maximum values Xmax depend on input values which are given in Tables 5 and 6, one concludes that from two 
Tables, the maximum value of downwind distance depends on stability classes and the effective height of point source, one finds 
that from the two Tables, there is well agreement between maximum downwind distance in actual and estimated values. 


Table 5 accuracy for the calculation of Xmax by using power law function 


H=5 H=100 H=250 
- H=5 H=5 . H=100 H=100 ; H=250 H=250 ; 
Stability ; Relative : Relative . Relative error 
Actual Estimated Actual Estimated Actual Estimated 
Class error (%) error (%) (%) 
Xmax(M) Xmax(M) Xmax(M) Xmax(M) Xmax(M) Xmax(M) 
A 11 11 0 300 287 4.33 1000 786 214 
B oi a 0 300 287 4.33 1000 786 214 
Cc 14 16 14.3 500 513 2.6 2000 1489 25.55 
D 45 35 22.22 1400 1637 16.92 5000 5301 6.02 
E 400 311 22.25 15000 21174 41.16 13000 76966 492 


Table 6 accuracy for the calculation of Xmax by a standard method 


H=5 H=100 
= H=5 H=5 . H=100 H=100 ; H=250 H=250 H=250 
Stability . Relative p Relative : ; 
Actual Estimated Actual Estimated Actual Estimated Relative 
Class error (%) error (%) 
Xmax (m) Xmax(M) Xmax(M) Xmax(M) Xmax(M) Xmax(M) error (%) 
A 0.482 0.476463 1.15 0.483 0.482756 0.05 0.483 0.482957 0.005 
B 3.50 3.34225 4.5 3.64 3.65818 0.5 3.66 3.66913 0.25 
C -2.90 -2.94464 1.54 -2.78 -2.78567 0.2 -2.78 -2.78093 0.03 
D -1.51 -1.55376 2.9 -1.52 -1.52362 0.24 -1.52 -1.52269 0.17 
E -1.71 -1.75254 14.2 -1.71 -1.725 0.87 -1.71 -1.7247 0.85 


In power law and standard method, one finds that the percentage of the error between estimated and actual maximum 
downwind distance in the range of 0 to 492 and 0.005 to 14.2 respectively. This work shows that when we used the method of the 
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power law and standard method are better than Brigg’s model method with respect to the previous relative error (Ronbanchob 


Apiratikul 2015). Then one can deduce that the best values for the estimation Xmax be observed at the lower Xmax value. 


=—@— actual Xmax at == actual Xmax 
5m at 100m 


=f estimated == estimated 
Xmax at 5m Xmax at 100m 


=—@— actual Xmax at 250m 


=—f— estimated Xmax at 250m 


Figure 1 The variation of the maximum estimated and actual downwind distance (m) using power law 


== actual Xmax at 250m 


== estimated Xmax at 250m 


==@ actual Xmax at == actual Xmax at 
5m 100m 

=fl— estimated == estimated 
Xmax at 5m Xmax at100m 


Figure 2 The variation of the maximum estimated and actual downwind distance (m) using the standard method. 
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One finds that the two figures show that the variation of the maximum estimated and actual downwind distance in power law 


and standard methods. From these figures, we conclude that the maximum estimated and actual downwind distances using 
standard method are better than using power law method but the power law method is a good agreement in all stabilities except 
stable condition. 


Table 7 shows Accuracy of the calculation of Cmax by using power law method 


ze H=5 H=5 H=100 H=100 H=100 Haae6 H=250 H=250 
Stability 7 Estimated Relative | Actual Estimated | Relative 7 Estimated | Relative 
Actualcmax1 Actual Cmax 
Class Cmax error(%) | Cmax Cmax error(%) Cmax error(%) 
0-*(g/m*) 10“(g/m”) 
10° (g/m’) 10“(g/m*)_| 10“(g/m?) 10°“(g/m*) 
A 96 96 0 0.24 0.24 0 0.035 0.035 0 
B 96 96 0 0.24 0.24 0 0.035 0.035 0 
C 84 85.8 2.14 0.21 0.21 0 0.0307 0.03 2.28 
D 60.3 64.37 6.74 0.1608 0.1609 0.06 0.0256 0.0256 
E 17 18.12 6.58 0.039 0.0453 15.38 9.21*107 9.21*107 
Table 8 shows the accuracy of the calculation of Cmax standard method 
H=100 H=5 H=5 H=100 H=100 H=250 H=250 
33 ; : H=100 . . H=250 ; : 
Stability Actual Estimated Relative Estimated Relative Estimated Relative 
Actual Cmax Actual Cmax 
Class Cmax10° Cmax10° error(%) Cmax error (%) Cmax error (%) 
10*(g/m?) 10“(g/m*) 
*(g/m*) *(g/m*) 10°“(g/m*) 10“(g/m*) 
A 1545.50 - 2950 290 159.253 42.783 73.13 1.56611 0.35898 77 
B 1212.28 1249.9 3.1 38.9136 41.563 6.8 23.627 -26.807 213 
C 368.313 324.175 11.98 16.3842 7.3607 55.07 1.27759 0.49106 61.5 
D 1625.53 -1804.4 211 68.5025 -8.1635 111.9 45.7823 0.000718 99.9 
E 1204.35 -1603.9 233 0.010702 -0.000968 | 109 2.7*10°9 -5.25 1.94*10%9 


Tables 7 and 8 show that the estimated of Cmax which depend on the maximum downwind distance, effective stack height, and 
stability classes. One finds that the errors are in the range of 0 to 15.38 and 3.1 to 290 for power law function and standard method 


respectively. 


4. CONCLUSION 


In this study, the dispersion parameter oy ando, can be calculated by using two methods the power law and standard methods at 
three effective height at 5, 100 and 250m respectively to calculate the maximum downwind distance and ground level concentration 
for the emission of air pollutants at point source based on Gaussian model. One finds that the percentage of the error between 
estimated and actual maximum downwind distance in the range of 0 to 492% and 0.005 to 14.2% respectively and the error in the 
previous work (Ronbanchob Apiratikul 2015) in the range of 0 to 2713.3%. This work shows that when we used the method of power 
law and standard method are better than Brigg’s model method with respect to the presented error. Then one can deduce that the 
best values for the estimation Xmax be observed at the lowest Xmax value. 

For the maximum concentration, one finds that the errors are in the range of 0 to 15.38% and 3.1 to 290% in power law and 
standard methods respectively. Also from the figures, we conclude that the maximum estimated and actual downwind distances 
using standard method are better than in power law method although the power law method is a good agreement in all stabilities 


except stable condition. 
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